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The chiral magnetization properties of cold and hot vacua are studied using quenched simulations in lattice Yang-Mills theory. In 
weak external magnetic fields the magnetization is proportional to the first power of the magnetic field. We evaluate numerically 
the coefficient of the proportionality (the chiral susceptibility) using near- zero eigenmodes of overlap fermions. We found that the 
product of the chiral susceptibility and the chiral condensate equals to 46(3) MeV. This value is very close to the phenomenolog- 
ical value of 50 MeV. In strong fields the magnetization is a nonlinear function of the applied magnetic field. We find that the 
"nonlinear features of the magnetization are well described by an inverse tangent function. The magnetization is weakly sensitive to 
temperature in the confinement phase. 
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'~J^ 1. Introduction 

Oh 

^ Quark is an electrically charged spin-1/2 particle with a mag- 
1^ I netic moment. Following an analogy with electrodynamics jlj], 
one can conclude that there exist at least two opposite types 
'of the quark magnetic back-reaction of the QCD vacuum on 
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the external magnetic field. The paramagnetic effect enhances 



QQ magnetic field in the vacuum due to polarization of the mag- 
netic moments of virtual quarks by the external magnetic field. 
On the contrary, the diamagnetic effect weakens the external 
\^ field due to the (quantized) transverse orbital motion of the vir- 
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tual quarks. 

Qualitatively, there is a distant similarity of the QCD vac- 
uum to a gas of electrically charged magnetic dipoles in elec- 
trodynamics f^. However, there also is an important difference 
between QCD and electrodynamics: in the physically relevant 
limit of two massless quarks, QCD is classically conformal the- 
ory in which mass scale appears due to nonperturbative quan- 
tum effects. The response of the vacuum on the external mag- 
netic fields is characterized by dimensionfull quantities and this 
feature makes the problem essentially nonperturbative. 

In our paper we study, for the first time, the magnetization 
properties of the non-Abelian vacuum using numerical simula- 
tions in lattice gauge theory. The problem is not limited by aca- 
demic interest, because very strong Abelian magnetic fields can 
be created, for example, in heavy ion collisions |3]. Such fields 
may substantially modify the phase diagram of QCD, changing 
even the order of the phase transition from the hadron phase to 
the quark-gluon plasma regime ^ . In lattice QCD the external 
field methods were already used to evaluate the electric polar- 
izability of neutral mesons and baryons [5J. The dependence of 
the quark condensate on the strength of the uniform magnetic 



field was also calculated recently in lattice gauge theory (|6|]. 

We concentrate on the (para)magnetic response of the vac- 
uum which appears due to polarization of the spins (magnetic 
moments) of the virtual quarks and antiquarks in the external 
electromagnetic field. A natural quantitative measure of the 
spin polarization in the vacuum is given by the expectation 
value 



where 

1 



-.X{F)mqFap 
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is the relativistic spin operator, Fap is the external electromag- 
netic field strength tensor and q is the electric charge of the 
quark For simplicity, we omit flavor indices in ([TJ and con- 
sider it for one quark flavor 

The quantity ([T]i was first introduced by lofFe and Smilga 
in Ref. fl] in order to analyze the nucleon magnetic moments 
which are related to phenomenologically interesting radiative 
transitions. Later the value of the magnetic susceptibility was 
estimated using various analytical approaches Jsl |^ [13, 111- 
The value of the magnetic susceptibility can be measured in ex- 
periments on lepton pair photoproduction through the chiral- 
odd coupling of a photon with quarks Iil21 . and in radiative 
heavy meson decays lll3ll . 

The right hand side of Eq. ([TJ is proportional to the elec- 
tromagnetic field strength tensor F^y - df,ay - dvUft due to 
the Lorenz covariance. The quark electric charge q appears 
in Eq. ([TJ since the electromagnetic field interacts with the 
quark field only in the combination qo/j. Another proportion- 
ality factor in the right hand side of Eq. ([TJ is the chiral con- 
densate ("P^P) (evaluated at the external electromagnetic field 
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F). This factor allows us to disentangle nonlinear effects of the 
enhancement of the chiral condensate in the external magnetic 
field lEIl from the effects of the quark's spin polarization. 

In a leading order the magnetization of the QCD vacuum in 
weak magnetic fields should be a linear function of the field 
strength. Using perturbative QCD and the Schwinger proper 
time formalism one can show that the magnetization due to 
a strong magnetic field at one-loop order is proportional to 
filogfiill. 

The strength of the vacuum polarization is characterized by 
a chiral magnetic susceptibility xiF) in Eq. ([TJ. In our discus- 
sion we treat the vacuum magnetization and the quark's spin 
polarization on equal footing skipping the ^-factor which re- 
lates these quantities and characterizes the gyromagnetic ratio 
of the quarks. 

The (chiral) magnetization of the QCD vacuum in the exter- 
nal magnetic field B = F12 - -F21 can be described by the 
dimensionless quantity 



so that 



(3) 



(4) 



for the other Lorentz components the polarization ([U is zero. 
The quantity (|3]l is of central interest in our paper. 

In Section |2] we derive an analytical formula which exactly 
relates the spinor structure of the low-lying Dirac eigenmodes 
to the magnetization ((S). This formula is used to evaluate the 
magnetization numerically in Section |3] We discuss in details 
both the linear magnetization in a weak field and nonlinear fea- 
tures of this quantity in stronger fields. Our conclusions are 
summarized in the last Section. 



2. A magnetization analog of Banks-Caslier relation 

In order to calculate the polarization properties of the QCD 
vacuum we derive the analytical formula which relates the mag- 
netization to the spin structure of the low-lying quark eigen- 
modes. This formula is an analogue of the Banks-Casher rela- 
tion lET 



<<P>P) = - lim 



np{A) 
V 



(5) 



which relates the chiral condensate (*I'*I') to the expectation 
value of the spectral density of the Dirac eigenmodes p{A) (to 
be defined below). In Appendix A we derive the Banks-Casher 
relation ^ in order to illustrate its relation to our analytical 
result (T9[ . 

The Euclidean partition function of QCD is given by the inte- 
gral over the gluon fields A^, yU - 1 , . . . , 4 and a - 1,. . .,N^-l, 
and over the quark Dirac fields 'V j, f - I, . . .Nf, 



Zqcd = JdaJd^'J D'I'e-^™(^>-^^('*'*''^> 
= J DA det[^(A) + M] e''^™^^^ , 



(6) 



where S ym is the Yang-Mills action. The fermion action is 

5f(A, >!',?') = Jd\^'f(x)[P(A) + M]ff,'i'f,(x), (7) 

where Ai is the Nf x Nf mass matrix in the flavor space. 

For the sake of simplicity we consider below one fermion 
species (Nf = 1) with the mass m. Then Ai = m and the 
fermion determinant in Eq. ^ is 



F(m) = det[^(A) + m] = Y\iiMA) + m) , 



(8) 



where - A^iA) is the eigenvalue of the massless Dirac oper- 
ator P = y/jOf, in the background of the Euclidean gauge field 
configuration A. The spectrum of this operator is defined by the 
equation 



(9) 



where ij/k — i//ii(x;A) is the corresponding eigenfunction (below 
we omit the argument "A" in i/zic and Ak for the sake of simplic- 
ity). Due to the anticommutation property, jsP + Pys - 0, any 
nonzero eigenvalue, A/t + 0, comes in a pair with its opposite, 
A-i - -All, corresponding to the eigenfunction tfr^k - js^k- 

The eigenfunctions (/r^ form a basis in the spinor space: they 
are orthonormalized and complete (we always omit spinor in- 
dices). 



^ 4fk(x)tf/l(x') = 6(x - x') . 



(10) 
(11) 



In the thermodynamic limit the expectation value of the mag- 
netization O) can be expressed via the nonzero eigenmodes, 
ip^k(x), of the massless Dirac operator (|9]l. 



-<>I'2„^>P) = <Tr[2,^>P(x)^(x)]) 

l//l(x) Ya/Bl^kix) 



= (Tr[2„^D(x,x)]) = <^ 

* 

iAk + m 



iAk + m 



At>0 



(12) 



<Zj — ^' 

,ii->0 



-iAk + ni 



where the trace is taken over the spinor indices and the fermion 
propagator D{x, x') is defined by Eq. ( l43T l. In the first raw of 
Eq. (fT2l l the notation (...) means the average over the gluon, A^, 
and fermion, 'V, fields. After the integration over the fermion 
fields is done, the brackets (. . .) mean the average over gluon 
fields with the weight F{m)e^^™'^^\ where F{m) is the deter- 
minant of the fermion operator ([8]l. We omit the explicit de- 
pendence of the eigenfunctions Ak and the eigenvalues i/r^ on 
the gluon field A^ and the external electromagnetic field B. In 
the first raw of Eq. (fT2l) the minus sign appears due to anticom- 
mutative nature of the fermionic fields. Following the logic of 
the derivation of the Banks-Casher formula, we ignore in sums 



over eigenvalues the exact zero modes, because the zero modes 
are inessential in the thermodynamic limit. 

The spin operator (|2| commutes with the 75 matrix, 

ys'^afi - '^afiys ^ , (13) 

so that ysSff/sys = because - 1- Then Eq. ( fT2] l gives us 
the following expression for the magnetization: 



{"V^ap^) = 2m{2_^ -5 ) . 



(14) 



Similarly to the derivation of the Banks-Casher relation (Ap- 
pendix A) we take the limit m — > and get 
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(15) 



where v(A) is the spectral density of the Dirac eigenvalues, 

p(A) = <v(^)) , v{A) = J] <J(^ - h) . (16) 

k 

Then we take the average of this expression over the whole 
space-time and take the integral over A 

nv{A) 



{^Y,p^) = lim<- ^ 



(17) 



Note, that if in Eq. ( flTI i we take the unit operator H instead of 
the spin operator Y.ap then we immediately recover the Banks- 
Casher formula Q due to the normalization condition (fTOl i. 

In Appendix B we give some arguments in favor of the valid- 
ity of the factorization (fTTl ). Then, using (|5]l, we get: 



(18) 



Our sketch of the proof of factorization property (fTsT i in Ap- 
pendix B is valid in infinite volume. In the next Section we 
check numerically the factorization in our finite volumes as 
well. 

Having compared Eq. ( fTSl l with the definition for the mag- 
netic susceptibility Xw in Eq- O, we get 

XiqF) qF„p = - Um< f d\ ifrlix; F) i/zAx; F)). (19) 

i->0 J 

Here (/'^(x; F) is the eigenmode of the Dirac operator in the ex- 
ternal (magnetic) background field B = Fi2. Next, we use the 
definition (O of the magnetization jj to rewrite Eq. (fT9] l as fol- 
lows 



= - iim< r 



d\ilj\{x-B)Zn^lfA(x;B)). 



(20) 



This is our final analytical expression which we use for the eval- 
uation of the magnetization in our numerical simulations. 

Eq. ( fTSl l demonstrates the apparent factorization of the chiral 
condensate in line with the original definition ([TJ of Ref. 
Due to the factorization (fTSl l the chiral condensate does not 
enter explicitly our final formulas (fT9] l and ( |20] |. and there- 
fore one can hope that various ambiguities (related to the def- 
inition and/or logarithmic divergence of the condensate in the 
quenched limit) does not enter our definition of the magnetiza- 
tion. 



Table 1: Parameters of simulations. 
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3. Magnetization from first principles 

3.1. Numerical simulations 

We simulate lattice SU(2) Yang-Mills theory with tadpole- 
improved Symanzik action for the gauge fields III8II . This 
improvement provides us with smoother gauge configurations 
compared to the usual Wilson action. The fermionic eigen- 
modes are calculated using the overlap Dirac operator for 
quarks in the fundamental representation [19]. One of the most 
important advantages of the overlap Dirac operator is its ex- 
plicit chiral symmetry at all lattice spacings. Basically, we used 
the same technique as was implemented in Ref. [,20i1 . 

In our quenched simulations we calculated the magnetiza- 
tion of the t/-quark condensate which has the smallest (absolute 
value of) the electric charge q -\q\ - elli. In continuum nota- 
tions, the uniform magnetic field in the third spatial direction is 
introduced into the Dirac operator by shifting the non-Abelian 
vector potential by the singlet Abelian potential a^: 



(21) 



BI2[x25^ 



I0fi2 



In order to adopt the field to a finite volume with peri- 
odic boundary conditions we have introduced an additional x- 
dependent boundary twist for fermions following Ref. 02111 . In 
a finite spatial volume with periodic boundary conditions the 
total magnetic flux trough any two-dimensional face of the 
lattice cube should be quantized |22]. This condition leads to 
the quantization of the uniform magnetic field: 



Ink 



keZ. 



(22) 



The physical strength of the magnetic field is a periodic func- 
tion of the flux number k with the period L^. The maximal 
strength is reached atk- 1? 12. 

The parameters of our simulations are given in Table I. For 
r = we used three lattice volumes and two values of the lat- 
tice spacing a in order to check the systematic errors due to 
finite volume and finite lattice spacing. In our simulations we 
used 20 gauge field configurations for each set of parameters 
given in Table I. We make simulations at zero temperature and 
at r = 0.827"^. The critical temperature m SU (2) gauge the- 
ory is Tc = 313.(3) MeV Q. We evaluate the Umit ^ ^ in 
( l20b by averaging this expression over low-lying nonzero Dirac 
eigenmodes with eigenvalues in the interval [0 . . . 50 MeV]. In 
order to evaluate the magnetization we use the asymptotic for- 
mula ( l20l i which is based on the factorization property ( fTSl l. 
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This approach is valid in the thermodynamic and chiral Hmits, 
taken simuhaneously. 

We show the magnetization at zero temperature as the func- 
tion of the external magnetic field in Figure [T] The values of 
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Figure 1 : The magnetization ,u as a function of the strength of the external 
magnetic field qB at zero temperature. The lattice spacings and spatial volumes 
are shown for all data sets. The fits are discussed in the text. 

the magnetization obtained at different spatial volumes and dif- 
ferent lattice spacings are very close to each other The rela- 
tive discrepancies are much smaller than, e.g., for the values 
of the chiral condensate |6]. This fact indicates finite-volume 
and finite-spacing dependencies of bilinear fermionic operators 
cancel in (|4]i with a good precision. 

In Figure|2]we have also checked that the factorization prop- 
erty (fTSl l is established very well for all checked strengths of the 
magnetic field. The values of the magnetization, evaluated with 
the help of the original (fT4l) and factorized ( l20l i definitions agree 
with each other within error bars. Unfortunately, the nonfactor- 
ized definition (fT4l i gives us much larger error bars compared to 
the factorized definition (|20] | at the same statistics. Therefore 
below we use the factorized definition only. 
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Figure 2: The check of the factorization property for the magnetization ft. The 
magnetization is calculated using both the oiiginal 1141 and the factorized )20) 
definitions (shown as the empty squares and the full circles, respectively). 

The behavior of the magnetization is consistent with general 
expectations: at low magnetic fields the magnetization is linear. 



This fact indicates the existence of a nonzero susceptibility at 
vanishingly small external magnetic field confirming the pres- 
ence of the paramagnetic contribution due to the quark's spin. 
At high magnetic fields the quarks ensembles should be fully 
polarized so that the magnetization should come to a saturation 
regime. Mathematically, in Eq. ( |20] | at strong magnetic fields 
the lowest-lying eigenfunctions with m - + \ eigenvalue of the 
spin projection operator E12 become dominated over the eigen- 
functions with the m - -\ eigenvalue, thus leading to the full 
polarization of the eigenmodes. In our units the saturation con- 
dition is to be as follows: 



lim niqB) 



-1 . 



(23) 



Figure [T] supports this observation. 

In Figure [3] we compare our zero-temperature data with 
the finite temperature magnetization obtained alT - Q.%2Tc- 
Visually, the effect of the temperature on magnetization is rather 




qB, GeV-^ 

Figure 3: The magnetization /j as a function of the strength of the magnetic 
field, qBalT = a (14* lattice) and T = 0.827",. (16^ x 6 lattice). Best arctan- 
based fits )30) with two fitting parameters are shown. 



small. More detailed analysis of T = and T > Q data will be 
done in the next subsection. 



3.2. Analysis of magnetization 

Let us compare our results with well known expressions for 
simple paramagnetic systems. The specific (i.e., per atom) mag- 
netization of a classical ideal paramagnetic gas is given by ||2[] 



cotha(B)- 



a{B) 



^lB 



(24) 
(25) 



where yu is the magnetic moment of a single molecule, T is tem- 
perature and ks is the Boltzmann constant. The function in the 
brackets in Eq. ( l24l i is often called the Langevin function. 

The quantum analogue of the Langevin function ( l24l ) for a 
spin 1/2-particle is described by the Brillouin behavior Il24ll : 



^quant(^) = ^ [2cOth2Q'(B) - COth Q'(B)] 



(26) 



The atoms are supposed to be electrically neutral otherwise 
the paramagnetic magnetization - either classical (l24l l or quan- 
tum ( l26b - should be supplemented by a diamagnetic contri- 
bution due to Landau quantization (the motion of the atoms in 
transverse direction to the magnetic field is to be constrained to 
the Landau levels). In an electron gas the diamagnetic contri- 
bution provides an essential part of the total magnetization HI]. 
However, we do not consider the diamagnetic contribution of 
the quark's motion because of the paramagnetic nature of the 
condensate ([l]). 

A similar behavior of the magnetization (l24l i and ( l26t is also 
provided by inverse tangent function which we add here for 
completeness: 



^'rig(B) ^ ^arctana(B). 



(27) 



At weak magnetic fields all functions ( |24] |. (|26] | and (l27T l are 
linear in the magnetic field. At strong magnetic field all expres- 
sions for the magnetization are consistent with the saturation 
property (|23] |. 

In the examples (l24l l. (|26] | and (|27] | the temperature T enters 
the magnetization via Eq. dZSl l. The temperature plays a role 
of the disorder ingredient which diminishes the magnetization 
of the system. In our simulations the disorder factor is coming 
- even at the zero-temperature case - from the random fluctu- 
ating background of the non-Abelian gauge fields, making the 
temperature T to be an eff'ective parameter in all considered ex- 
amples. 

We fit our data by three functions inspired by Eqs. dZTl l. (l24l l 
and ( 



..quant, „x 



- fio 



coth 



2 coth 



TTXoqB 



3xoqB 



coth 



xoqB 



2/Zoo 

arc tan 

71 Zfl. 



(28) 
(29) 
(30) 



These fitting functions share the following properties. All these 
fitting functions are the functions of the two parameters: the 
zero-field susceptibility xo and the strong-field saturation con- 
stant, yUoo. Indeed, in weak and strong field limits one gets, re- 
spectively 



fi^iB) ^ XO ■ qB + 



as <7Z? — > ; 



and 12511 

lim fi'f^fiqB) = -fi^ 



(31) 



(32) 



Here the index / stands for the type of the fitting function: 
"class", "quant" or "trig". 

According to Eq.(l23]l we expect that 



1 



(33) 



Thus we fit the data using the functions ( l28T l. ( |29] | and dSOl l with 
Xo and fioo being the two fitting parameters. We also perform an 
additional set of one-parameter fits with fixed /loo - I- Finally, 



we also fit the weak-field behavior of the magnetization by the 
linear function (ISTT i. 

The fitting results are presented in Tables |2] and [3] for zero 
and non-zero temperature cases, respectively. 

Table 2: Fitting results at T = 0. 
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Table 3: Fitting results at T = 0.821 
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-X(h GeV ^ yUco 


x'IA.o.t 




one parameter fits 




linear 


1.46(6) 


1.8 


arctan 


1.46(2) 1 


1.05 


Langevin 


1.36(4) 1 


3.8 


Brillouin 


1.19(6) 1 


17.6 




two parameter fits 




arctan 


1.53(3) 0.94(2) 


0.57 


Langevin 


1.49(4) 0.87(2) 


0.71 


Brillouin 


1.47(4) 0.71(2) 


1.1 



At r = all fits have quite large values of x^ lA.o.f., Table|2l 
The lowest value of this important quantity is reached for the 
arctan fitting function ( l30l l with two fitting parameters. Notice 
that both for the arctan-based ( l30b function and for the classical 
Langevin ( |28] | fitting function the asymptotic values of the po- 
larization yUoo are very close to the theoretical expectation (l33b . 

At r = the quantum (Brillouin) function ( |29] l does not 
work at all: it has the very large value of x^ l^-O-i- and the 
strong-field limit is also inconsistent with our expectation (l33T l. 
The fitting functions - corresponding to the one-parameter fits 
- are shown in Figure [T]by lines. 

We point out that the linear fit at the weak field limit does not 
agree with most other fitting functions. The linear fitting is done 
for the relatively weak magnetic fields, < ^qH < 500 MeV 
(shown as a straight line in Figure[T]l. The fields from this inter- 
val, however, are of the order or even higher than typical QCD 
scale Aqcd ~ 200 MeV. Therefore, the nonlinear eff'ects may 
affect the determination of the zero-field susceptibility if the fit- 
ting is done by the linear function. Below we quote the result 
for the magnetic susceptibility using the arctan-based ( |30] | func- 
tion. The corresponding value is presented in the bold font in 
Table |2] 

In a finite volume the physical magnetic field (l22l i is a pe- 
riodic function of the number of elementary magnetic fluxes k 
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going through any face of the lattice volume f2^. Since the 
periodic finite-volume behavior at very strong fields is not re- 
flected in the form of flie fitting functions (l28l l. ( |29l ) and dSOl l the 
asymptotic magnetization /Uoo of two-parameter fits may deviate 
from the expected high-field limit (l33T l. However, the artifacts 
related to the finite volume of the lattice are small since for the 
best fits by the functions (l28T l and dSOl l the deviations from the 
limit ( l33T l are small. We also note that the logarithmic effects 
- predicted in Ref. lfl4ll - cannot be reliably determined from 
our data because of the coarse grid of the data points at strong 
fields. 

At r = 0.82rc the fits have more reasonable values of 
X^/d.o.f., Table[3] The error bars of the best fit parameters are 
larger so that the zero-field susceptibility agrees within error 
bars for almost all the fitting functions except for the Brillouin 
function (|29] ). The best fit - in terms of both the quality of the 
fit and the value of the asymptotic polarization jjca, Eq. ( l33T l - is 
the arctan-based (|30] | function. The accepted value is given in 
the bold font in Table |3] Finally, the linear fitting is done in the 
range < y/qH < 0.425MeV. 

The best two parameter arctan-fits for both values of temper- 
ature are shown in Figure [3] by lines along with the data. The 
magnetic susceptibility turns out to be insensitive with respect 
to the variation of the temperature: 



-1.547(6) GeV" 
-1.53(3) GeV" 



r = 

T = 0.827^ 



(34) 



Our lattice spacings, a - A^^, correspond to the scales 
Auv(r = 0) ~ 2GeV and Auv(r = 0.82 T,.) ~ 1.5 GeV, re- 
spectively. Other estimations of the chiral susceptibility were 
done in 

Theoretically, the value of the magnetic susceptibility can be 
parameterized in the form li9J 



(35) 



were is a dimensionless parameter and f„ = 130.7 Me V is 
the pion decay constant for Nc - 3. 

In the notations of Eq. ( |35] | the result of our calculation (l34l l 
at r = corresponds to 



c (T = 0) = 1.043(4). 



(36) 



The operator product expansion combined with the pion dom- 
inance idea gives us the value = 2 The corresponding 
theoretical prediction for S U (2) gauge theory is 



;^.(^ = 0) = -2.97 GeV 



-2 



(for A?,. = 2). 



(37) 



The holographic description of QCD gives a slightly higher 



value for the su scep tibility, Cj, = 2.15 lllOll . Both the results 
of Refs. 10] and jw) agree well with the original QCD sum rule 
fit made by lofife and Smilga [7]. 

In the instanton liquid model of the QCD vacuum the mag- 
netic susceptibility was first calculated in Ref. fill]. After a 
proper rescaling we obtain the value - I -24 corresponding 
to the instanton model. 



Thus, our zero-temperature result ( |34] |. ( l36l l is by 25% 
smaller than the value of the magnetic susceptibility obtained 
in the instanton vacuum Ol ill , and is by a factor of two smaller 
than the value obtained by traditional field theoretic and mod- 
ern holographic approaches. These discrepancies are not un- 
expected since we used the quenched lattice study in which all 
vacuum quark loops are ignored, and the anomalous dimension 
of the chiral susceptibility was not taken into account. More- 
over, our calculations are performed in the S U (2) gauge theory 
in which the number of colors is reduced in comparison with 
the real QCD. 

An experimentally relevant and phenomenologically inter- 
esting quantity is given by the product of the chiral suscep- 
tibility X and the chiral condensate (*I'*I') 112|]. Using our 
zero-temperature result for the chiral susceptibility (|34] | and the 
value for the chiral condensate obtained in other quenched stud- 
ies [6,2Qh = [310(6)MeV]\ one gets 



- X = 46(3) MeV [quenched limit] . 



(38) 



This result is surprisingly close to the estimation based on the 
QCD sum rules techniques, which gives for ( l38T l the number of 
tiie order of50 MeV in. 



4. Conclusions 

We have evaluated for the first time the magnetic susceptibil- 
ity of the chiral condensate using the first-principle methods of 
lattice S U{2) gauge theory in the quenched limit. To this end 
we have derived formula (fT9] l which relates the chiral magneti- 
zation of the QCD vacuum to the low-lying chiral eigenmodes 
of the Dirac operator in a manner of the Banks-Casher relation 
(the exact zero modes do not contribute to the magnetization in 
the thermodynamic limit). In order to derive Eq. (fT9] l we used 
the factorization property ( fTSl l. which was verified numerically. 
We calculated these eigenmodes using the overlap fermion op- 
erator in the background of gluon fields generated with the help 
of the tadpole-improved Symanzik action. 

We found that at weak magnetic field the magnetization of 
the QCD vacuum is a linear function of the field strength. The 
associated chiral magnetic susceptibility is almost independent 
on temperature ( l34l i up to T - 0.82 T^.. The value of the mag- 
netization (|34] | is smaller compared to the existing analytical 
estimates (l37l i for S U{2) gauge theory. We attribute the reason 
for the difference to the quenching effects. 

The nonlinear features of the magnetization are very well de- 
scribed by an inverse tangent function ( l30t of the applied mag- 
netic field. Figure [3] This parametrization of the magnetization 
works both at zero and non-zero temperatures in the confine- 
ment phase. 
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A. Banks-Casher relation: chiral condensate via eigen- 
modes 

In Section |2] the relation of the magnetization to the Dirac 
eigenmodes was shown to be given by Eq. ( fT9b . For the sake 
of completeness we present below a derivation of the Banks- 
Casher relation UtIi which is very similar to the derivation of 
Eq. (O. 

A differentiation of the partition function Q with respect to 
the mass m in the thermodynamic limit provides us with the 
chiral condensate. 



V dm 



In 



Ai,>Q k 



(39) 



where we have used Eq. (O as well as the pairwise appearance 
of the mutually opposite eigenvalues. 



1 dF{m) 
F(m) dm 



Zl 2m 
iAk + m f-^ + m 



(40) 



At>0 k 



The zero mode(s), corresponding to Aq - 0, are not counted in 
Eq. (|40] | because they give vanishing contribution in the ther- 
modynamic limit, and thus are irrelevant. 

In the chiral limit, m — > 0, Eq. (l39T l provides us with the 
celebrated Banks-Casher formula |17], Eq. (|5]l, which relates 
the chiral condensate (*I'M^) to the expectation value p(A) of the 
spectral density v(A) of the Dirac eigenvalues ( fTST i in the limit 
to zero virtuality, A ^ 0. In order to derive Eq. ^ we have 
used the relations 



lim 



1 



71 A^ +1 



6(A) and 



00 

/ 



dA6(A) 



(41) 



so that 



v-i 2m r 2m 
) -0 T= dAv(A)— 

00 

^271 J dAv(A)6(A)^7T]imv(A). 



(42) 



The fermionic propagator in the background of the gauge 
field Af^ can be expressed in terms of the eigenmodes (|9]l as 
follows: 



mxmx')h ^ D(x,x') = y ^.y^^' 

lAk + m 



(43) 



Due to the completeness condition (fTTl i the propagator (l43T l sat- 
isfies the equation 



[P(A) + m] D(x, x') = 6ix - x') . 



B. Factorization in Eq. dTSt 



(44) 



Below discuss a possible origin of the factorization in Eqs. 
([TtI i and (fTSl l. which was used to evaluate the chiral suscepti- 
bility in this article. Consider two local operators, Oiix) and 
02ix), which are functions of the lattice coordinate x. Then in 
infinite volume limit 



\^Yj^x{x)02{y))^{0,){02). 



(45) 



The proof of ( |45]) is trivial, since the quantum average in lat- 
tice calculations is equivalent to a sum over infinite number of 
gauge field configurations: 



(46) 



i^Oi(x)02(3')) 

X 



where V is the number of lattice points x, and the sum in ( l46b 
goes over gauge field configurations labeled by the index /. The 
gauge field configurations are generated with the Boltzmann 
probability density e"''™^'*'/^^- For infinite volume, V — > 00, 
the average for one gauge field configuration is equal to the 
quantum average. 



{Ox{x)). 



(47) 



Equation (|45] | follows from Eqs. (|46] | and (l47T l. 

If one takes 0\{x) - 02{x) - TrF^y(x), then in the contin- 
uum limit we get from Eq. (l45l l the well known factorization 
formula: 

limj^ i J d\ Tr F%(x) Tr F^Jy)"^ = <Tr F%f . (48) 

The factorization of the magnetization, Eqs. (fTTb and ( fTSl l, is 
very similar to the factorization i45[ . We only have to prove 
that the analogue of Eq. ( l47l i is valid for the bulk quantities 
lim^^(i7Tv(A)/V cindlimji^o J d'*xi/' j(x)2Qy3i/r^(x). One can nat- 
urally assume that these quantities have a smooth behavior to- 
wards the continuum limit, V — > 00. Then the average over one 
infinitely large configuration of the gauge fields is equal to the 
usual quantum average of the factorized expression, and then 
the factorization ( fTsT i is valid. 
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